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Abstract 

We study the hydrodynamic sound mode using gauge/gravity correspondence by 
examining a generic black brane background's response to perturbations. We assume 
that the background is generated by a single scalar field, and then generalize to the 
case of multiple scalar fields. The relevant differential equations obeyed by the gauge 
invariant variables are presented in both cases. Finally, we present an analytical solu- 
tion to these equations in a special case; this solution allows us to determine the speed 
of sound and bulk viscosity for certain special metrics. These results may be use- 
ful in determining sound mode transport coefficients in phenomenologically motivated 
holographic models of strongly coupled systems. 
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1 Introduction 



The plasma created in heavy ion collisions at RHIC appears to be both strongly coupled, 
and well described by hydrodynamics [TJ [2j [3]. Hydrodynamics is an effective theory which 
describes the behavior of a fluid on length and time scales which are much longer than any 
microscopic scale. The hydrodynamic stress-energy tensor is constructed from conserved 
quantities, and is built so that it respects equilibrium thermodynamics and the symmetries 
in the problem. One can then add plane wave type perturbations ~ ^{qz-wt)^ an( j uge 
conservation equation d p T^ u to analyze the normal modes of the system - see [U [5] for 
reviews. For the shear mode, the dispersion relation between the energy of the perturbation 
w and the momentum q is found to be 

w(q) shear = -ij—jsQ 2 , (!) 

while for the sound mode it is 

w(q) smin d = v s q - i—-^ + — )q 2 . (2) 

e + F \ p 2?7 I 

Here, p is the number of spatial dimensions, e and P denote the equilibrium pressure and 
energy density, respectively, v s is the speed of sound and r\ and ( are the shear viscosity and 
bulk viscosity. It is desirable to calculate these latter three quantities which are used in hy- 
drodynamic models of heavy ion collisions, but in the regime of strong coupling, perturbative 
calculations are unreliable. 

The AdS/CFT correspondence P,[TJ, E], (or more generally, the 'gauge/gravity duality' or 
'holography') has become an indispensable tool for describing strongly coupled systems, and 
has enjoyed arguably its greatest success calculating hydrodynamic transport coefficients. 
Within this framework, hydrodynamic dispersion relations are calculated by computing cor- 
relation functions of the stress-energy tensor, and then either examining the poles of such 
quantities, or by using Kubo relations [9j [TOj EH [121 EH E]. Alternatively, one can exam- 
ine behavior of a background under perturbations, and determine the dispersion relation 
by imposing appropriate boundary conditions [TBI, ITS] . The resulting dispersion relations 
from all of these methods coincide as demonstrated in [15]. After computing the dispersion 
relation, one can compare with the formulas ([1]), ([2]), and read off the appropriate transport 
coefficients. 

Another method which has been employed to access such quantities is the 'membrane 
paradigm', in which one computes the hydrodynamic properties of a black hole's 'stretched 
horizon'. These hydrodynamic quantities agree with the AdS/CFT calculations in all cases 
tested so far [TTJ HH [19], [201 EI]- The details of the membrane paradigm approach are not 
relevant for this paper, but we mention these works for historical reasons. In [IT], a formula 
for the shear viscosity to entropy density ratio rj/s was derived which is applicable to a 
wide variety of gravity duals. Application of this formula to known gravity duals always 
resulted in rj/s = 1/47T, and eventually led to the celebrated viscosity bound conjecture that 
v/ s ^ 1/47T for all physical substances. While the membrane paradigm approach was first 
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used to derive this general formula, it can also be derived by examining the quasi-normal 
mode spectrum as shown in [22]. Clearly, a similar general formula for the bulk viscosity 
would also be desirable in the hopes that other universal behavior might be discovered. Such 
a formula would also be useful when attempting to phenomenologically fit the results from 
the lattice [251123]. 

This work is motivated by the goal of developing such a formula, though such a formula 
will not be presented here. Instead, this work should be viewed as a first step towards this 
goal. Here, we consider sound mode fluctuations of a dual gravity theory supported by a 
single scalar field, and then generalize to the case of multiple scalar fields. We attempt to 
be as general as possible by not specifying the background profiles of the scalar field(s). 

The recent works J25J [26] also explore the sound mode in phenomenologically motivated 
single scalar models. The work presented here is complementary to these papers, though we 
allow for the possibility of multiple scalar fields, and also do not restrict ourselves to five 
dimensions. 

In section 1, we write the relevant equations that must be obeyed by the background 
fields. In section 2, we introduce the standard sound mode perturbations and explicitly 
compute the linearized Einstein equations. In section 3, we reduce these Einstein equations 
to two gauge invariant equations. In section 4, we illustrate how one uses these equations to 
solve for the hydrodynamical dispersion relation and thus compute the speed of sound and 
bulk viscosity. The gauge invariant equations are rather complicated, so we only consider 
a simple special case in this section. The results of this section are generalizations of the 
results of [16] . Finally, in section 5, we generalize our gauge invariant equations by including 
an arbitrary number of scalar fields. We conclude and mention some prospects for further 
investigation in section 6. 

Our general relativistic conventions are those of [27]. we use R^ u , G^ u and to denote 
the Ricci, Einstein, and stress-energy tensors respectively. We will also make use of the 
shorthand notation 

R° = F 



RZ = F„ 



X 



R r r = F r . (3) 

In Appendix [A[ one can find explicit forms of these functions in terms of the metric compo- 
nents. V M denotes the covariant derivative, and □ = V^V^. Throughout the paper, we also 
use the notation T>l to denote the logarithmic derivative, namely 

X'(r) 

Vl [X(r)) = ^ii (4) 

2 Background Fields 

We wish to examine hydrodynamic fluctuations on the following p + 2 dimensional gravita- 
tional background 

ds 2 = goo(r)dt 2 + g xx (r)dxjdx j + g rr (r)dr 2 (5) 
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where j = 1, 2...p, and r is the extra-dimensional coordinate. We will often make use of the 
definition 



f(r) = ^-g 00 (r)g**(r). (6) 

We assume the position of a horizon at r = tq, and that the behavior of the metric 
components near the horizon is 

WO ~ -7o(r - r ) + D(r - r ) 2 

ffrr(r) « -^-+0(1) 

r — r 

&r*(r) « fe(r ). (7) 

The quantities 70, 7 r and ^(ro) are independent of r. The Hawking temperature is given 
by 

T = ±M. (8) 
47r y 7 r 

Let us make the simplest assumption, that this metric is created by a single, minimally 
coupled scalar field. In other words, we assume the action is of the form 

S = ^Jq J d p+2 xV=g~ (R - \d,<t>d^ - U{4>)) . (9) 

Then, the background equations are 

G^y = SnGp^T^ (10) 
dU , x 

Here we ignore any subtleties regarding boundary terms that come from integration by parts. 
Such terms can be taken care of by adding additional boundary terms to the action, though 
we have not explicitly written such terms above. 

The energy-momentum tensor derived from the action is 

SttG p+2 T^ = \{d^d v <P-g^) (12) 
£<» = \d x <t>d x <t> + U{<l>) (13) 

Because the background metric only depends on the extra dimensional coordinate r, it is 
clear that the field <fi is also only a function of r. In the case of multiple scalar fields, it 
might be possible to have fields which depend on the other coordinates, provided that all 
such dependence cancels out in the combination that appears in T^ u . We will not consider 
such special cases, and will always assume that the scalar fields only depend on r. 
In this case, we can get a more explicit form of the potential: 

dUm) =nt = ^d r \^ 99 ^>(r)}. (14) 



d(t) 
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Throughout this work, the prime denotes derivative with respect to r. It will also be useful 
to express the second derivative of the potential as 



d 2 U(0(r)) 



-d r 



dU{(j>{r)) 



(15) 



dcj) 2 (j)'{r) 

One can write the field 0(r) in terms of the metric components by noting that 

g m G m - g rr G rr = F Q (r) - F r (r) (16) 



and 



8nG p+2 (g 00 T 00 - g rr T rr ) = \ [-£+ - (g rr <t>'{rf - C^r)) 



[IT) 
18) 



Thus, 
It 

seen by considering 



( j>>(r) 2 = 2g rr (F (r)-F r (r)). (19) 
It is also noteworthy that this background has the special property F = F x as can be 



g 00 G 00 -g xx G tt = -8nG p+2 (g 00 T 00 - g xx T u ) (20) 
F {r)-F x {r) = C^r) - C^r) = 0. (21) 
Because of this fact, and the general theorem given in [28], all backgrounds we consider have 

f]/s = 1/47T. 



3 Hydro dynamic Fluctuations 

We now introduce fluctuations of the fields on this background g^ u — > g^ u + 5g^ u and — > 
+ 5(fi, and assume the usual time dependence 

5 9ia ,(t,z,T) = e^-^h^r), (22) 
8(t>{t,z,r) = e i{qz ~ wt) 5(j)(r). (23) 

Here we use the coordinate z to denote one of the spatial coordinates: z = x p , and w and q 
are the energy and momentum of the perturbation. For the sound mode, in the gauge where 
h^ r = 0, the only non-zero fluctuations are [TTJ [16] 

h 00 (r) = goo(r)A(r), (24) 

1 

Y,hu(r) ee g xx (r)B(r), (25) 
P 1 i=i 

h zz (r) = 9xx(r)C{r), (26) 

h 0z (r) ee g 00 (r)D(r), (27) 

5<f>{r). (28) 

Turning on these perturbations, and expanding the background equations of motion to first 
order in the perturbation, we end up with a set of 8 differential equations for the perturba- 
tions. The equations are presented in full in Appendix IB1 
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4 Gauge Invariant Equations 

There is still residual gauge freedom in these equations under the infinitesimal diffeomor- 
phisms 



(29) 
(30) 



for any vector £ M = ^(r)e^ qz wt \ The following gauge invariant combinations given in [T6] 
transform only into themselves under such a diffeomorphism. (i.e. Zq — > and — > Z^.) 



(31) 



Z (r) = -/(r) 2 (q 2 A{r) + 2gu; J D(r)) + w 2 C(r) - q 2 ^f^- + w 2 B{i 

\ 9xx\ ) J 

<f)'(r) 



ZJr) 



5(f)(r) 



tB r) 



(32) 



By taking an appropriate combination of the equations (17211 79 p . one can arrive at two coupled 
second order equations for the gauge invariant combinations. The particular combination 
of the Einstein equations used, as well as additional details regarding this derivation, are 
presented in Appendix A combination which appears frequently in these equations is 



a(r) = q 2 (p- 1) 



V L [g xx (r)}J /W 2 ' 



(33) 



The following is the most compact form of the gauge invariant equations that we have 
been able to find. 



Qrr 2 M • 



-9 



--a 2 td r 



1-99 ry, 

a 2 / 4 
IZ^'f (ad,, 



+ Z (V L [f 2 }V L [f 2 a}-g rr (w 2 g 



,2„00 



2 -V 



pV L [g xx \ 



q 2 g xx 



-99' 



(34) 



-.d r 



-gg"Z'^ +-d r 



^l[9xx]\ \f" r [ f 



1 B ' 



{w 2 g m + q 2 g" 



9xx%4 



-z. 



(V <T + w 2 g 00 



(^L[g X x})\ 



f 



f 



2 AJ 



_(V L [gxx\y 



-99 



(35) 







Of course, one can write these equations only in terms of the metric components because of 
the relation (JT9"]1 . 

To determine the dispersion relation, one needs to solve the above equations perturba- 
tively in w, q, and then apply the incoming wave and Dirichlet boundary conditions [15]. We 
will illustrate this procedure in the following section. Some terms can be neglected when 
considering these equations to 0(w, q), but one should keep in mind that Z is 0(q 2 ), Z^ is 
0(1), and a is 0{q 2 ). 
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5 Example Application 



As a simple example of how to apply these equations, consider the special case of 0(r) 
C\ log [^(r)] where C\ is a constant. In this special case, the equation (135]) reduces to 



-fj 



^g rr Z^} - Z^g rr (q 2 g™ + w 2 g 00 ) = 0. 



This can be more easily seen by re-writing 



0. 



_(D L [g xx }) 
P 



T> L [g xx 

^9rrf 



2 V L 



'-99" (, 
^9 rr V L [g 



F> L [g xx ] 

^9rrf 
Vl[9xx\ 



(Fq - F x ) - f 2 d r 



(F - F x ) - pd r 



2g rr (F — F r 

. p(D l [9xx}) 2 . 
<P\r) 2 

p(T> L [g xx ]Y 



(36) 

(37) 

(38) 

(39) 

(40) 
(41) 



The standard way of applying the incoming wave boundary conditions is to make the 
ansatz [T5J [16] 



Z (r) = f(r)'^(Y (r) + qY 1 (r) + ...) 
Z*{r) = f(r)-^(Y <po (r) + qY^(r) + ...) 
w(q) = Wiq + w 2 q 2 + .-- 



(42) 
(43) 
(44) 



with the condition that all Y functions are regular at the horizon. Inserting this ansatz into 
f[3"6"j) . expanding the result in powers of q, and neglecting terms of 0(q 2 ) and higher, we find 
the following equation: 



^99^0 +qd r 



-99' 



YL - ^V L U)Y^ 



4>i 



2nT 



(45) 



Solving this order by order in q is now quite simple. The solution for Y"^ can be written in 
terms of an integral 

9rr(r') 



Ytpoir) = c + ci / dr 
Jro \ -9(r') 



(46) 



but this integral is logarithmically divergent at the horizon by ([7]). Thus, the assumption of 
regularity on the Y functions leads to Y^ — c o- Finally, this constant must be set to zero by 
the Dirichlet boundary condition at infinity. Plugging Y^q = into the next order equation, 
one also finds Y^i = 0. 

Next, one must solve the equation for Z using these boundary conditions with the 
knowledge that = 0. But before doing so, it is useful to pause and ask what type of 
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metrics these results will be applicable to. There are 3 unknown metric functions, g 00 , 9xx, 9rr, 
but we have two constraints on them, namely (1191) and (1211) . It can be shown using these 
constraints, and the assumption of the scalar field profile <fi ~ log(g xx ) (as at the beginning 
of this section), that we are considering metrics which satisfy Fq = F x and the following 
constraint: 

9oo(r) = -^-^{r) + a l9xx (rY^ +l . (47) 
a 2 -p 

Here, a , a%, a 2 are independent of r, and we have defined the coefficients as above for future 
convenience. Clearly, a\ can be determined in terms of ao and a 2 , by requiring g 0Q vanish at 
the horizon. 

Returning now to the equation for Z , going through the same steps of inserting the 
incoming wave ansatz, expanding in powers of q, solving order by order in q, and applying 
the boundary conditions leads to the dispersion relation 



w{q) = {-^ q - l ^r P q2 + <D{(f) - (48) 

The details of this derivation are presented in the Appendix [D] Comparing this dispersion 
relation with the expected hydrodynamic dispersion relation ([2]) yields the relations for the 
speed of sound and bulk viscosity: 



/ ao — a>2 



V p 

C 2(2a 2 -p + l) 



(49) 



Let us now explicitly check these results agree with known cases. Consider the 'Schwarzschild 
AdS Black Hole' metric in p+2 dimensions, 



d ° 2 - 15 



r i Tj^ dv^ 

f(r) 2 dt 2 + d Xj dx j + —. (50) 

r J{ r ) 



f(r) 2 = 1-(7) P+1 (51) 

This is equivalent to the metric (f47I) with the choices ao = (p+ 1)/2, a 2 = (p — 1)/2. Inserting 
these into (149]) yields 

v s = i/Vp 

C = (52) 

which are in perfect agreement with our expectations, since the metric is conformal, and is 
thus dual to a conformal field theory. 
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Next, consider the 'Dp-Brane metric' in the Einstein frame, which can be reduced to 

[HI EE] 



ds 



2 



l) ' [~f(r) W + dx^] + ( T T ) P (53) 



L) f(r) 



r \ 7 -P 



f(rY = ■ (54) 

This is the same as (jlTj) with the choices a = (7p — p 2 )/(9 — p) and a 2 = 2p/(9—p). Inserting 
these into (119]) gives 




c = 

7] 9 — J9 




(55) 

tj — p 

which is in agreement with the result of [16] . 

It has been conjectured [29] that the relationship 

(56) 

should hold for a strongly coupled plasma in p spatial dimensions. From our formulas, one 
can see that this relation is satisfied provided 

«o + «2 > P- (57) 

Both special cases considered above have ao + a 2 = p, and thus saturate the conjectured 
bound. At first sight it appears that it would be possible to violate the conjectured bound 
with an appropriate choice of metric, but whether or not such a gravity dual could be 
embedded into string theory, and whether or not such a model would provide a reasonable 
description of a strongly coupled plasma should perhaps be investigated in the future. 

This completes the example calculation for the single scalar field profile chosen. In the 
process of this example, we have generated equation (jl9l) . which is applicable to metrics 
which obey Fq = F x , and (HTj) . Equations (H9|) are generalizations of formulas given in [16] . 

This is only one (particularly simple) application of the gauge invariant equations in 
the previous section. In the next section we increase the generality of the gauge invariant 
equations by allowing for multiple scalar fields. 

6 Multiple Scalar Fields 

Consider now adding additional scalar fields, so that the total number is n. 

S = 16 J Q ^ J dP^x^Tg \R-\Y. d^kd^k - Ufa, <h...<f> n )J • (58) 
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There is now a background equation for each scalar field. As before, assume each field only 
depends on the radial coordinate r. 



□ c 



Adding these all together gives 

n n 1 

k=l k=l V 9 



-99 rr Ur) 



dU{(j)i,(j)2, •••0n) 



2-^ PtA. 



k=l 



In addition, the energy momentum tensor now becomes 

n 

k=l 

n 

£<f>k = \^d X <pkd X (t)k + U "(01, 4>2~4n) 
fc=l 

and the constraint equation ( fT9l) becomes 

n 

Y,<t>' k (r) 2 = 2g rr (F (r) - F r (r)) . 



(59) 
(60) 

(61) 
(62) 

(63) 



k=l 



It is fairly straightforward to generalize the gauge invariant equations to include more 
scalar fields. The linearized Einstein equations are modified due to the presence of other 
fields, and one needs to introduce an additional gauge invariant variable for each additional 
scalar field. We use the index i to denote a particular scalar field with i = l,2...n. 



Z<in{r) = 6(f>i(r) 



#(r) 



-B(r) 



£>l [9xx{r)\ 

One then gets additional terms in the differential equation for Zq. It can be written: 



(64) 



/rr 2 r4 



'-9 



a 2 f 4 d r 



'-99 



2/4 ^0 



£ |2Z # 0'J 2 [ad r 



+ Z (v L [f\V L [fa\ - g„ (w 2 g 00 + q 2 g 
1 (w 2 , 



+ 



a VP-" 



pT> L [g xx ] 



99 rr 4>' k 



(65) 



There are also n gauge invariant equations, one for each of the scalar fields. Each of these 
equations takes the form 

d 2 U 



9r 



n 



-9 



-gg rr Z^ - Z^g rr (w 2 g 00 + q 2 g xx ) - g rr £ Z, 



pT> L [g xx )a 
2 



E 

.k=l 



k=1 ^ dfcdfa 



af 



Or 



V L [g xx ] 



0, 



z <t>k<P' k \ 



rr if 

-99 <Pk 



+ pg X xV L 



v L [g* 



(w 2 g™ + g 2 g™ 



f 



0. 



(66) 
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At first sight, this equation appears slightly different than that of a single scalar field (135|) . 
There are two reasons for this. First, in the case of a single scalar field, the term d ^.g^ k 
can be written in terms of the metric only due to (|15|) ; one cannot do this if the potential 
depends on more than one variable. Secondly, in the case of a single scalar field we removed 
a term which vanishes by the equations of motion, but we have not done so here. The term 
in question appears in (1851) . 



7 Conclusion 

In this work we have presented a set of gauge invariant equations for sound mode perturba- 
tions on a generic black brane type background. The equations ( !65ll66l) are the main results 
of this paper. These equations can be used to determine the speed of sound and bulk vis- 
cosity for any metric which can be generated by a set of minimally coupled scalar fields. In 
order to determine the dispersion relation, one must solve these equations perturbatively in 
q, applying the incoming wave boundary condition at the horizon, and Dirichlet boundary 
condition at r = oo. 

The gauge invariant equations are quite complicated, and so far a general analytic solution 
has eluded us, though we did present a solution for a particular class of metrics. Metrics which 
obey (H7|) and F = F x have speed of sound and bulk viscosity given by (H9|) . These results 
are a generalization of the results of [16], and include both Dp-brane, and Schwarzschild AdS 
black hole metrics. 

In the future, we hope to report further on the possibility of an analytic solution of 
these equations. In the event such a solution is not possible, the equations can be solved 
numerically for a specified metric and set of scalar field profiles. One could also increase the 
generality of the gauge invariant equations by including other types of matter such as gauge 
fields. 

We believe that the equations derived here may be useful for phenomenologically based 
models such as those of [201 EI]- It would be interesting to see what these models (which were 
chiefly developed to match the meson mass spectrum) have to say about hydrodynamics. 



Acknowledgments 

I would like to thank Joe Kapusta for helpful comments. I would also like to thank the 
Physics Department at McGill University, and the organizers of the workshop "AdS/CFT, 
Condensed Matter, and QCD" for their hospitality. This work was supported by the US 
Department of Energy (DOE) under Grant. NO. DE-FG02-87ER40328, and by the Graduate 
School at the University of Minnesota under the Doctoral Dissertation Fellowship. 



11 



References 

[1] I. Aresene, et al. (BRAHMS Collaboration), Nucl. Phys. A757, 1 (2005); B.B. Back, 
et al. (PHOBOS Collaboration), Nucl. Phys. A757, 28 (2005); J. Adams, et al. (STAR 
Collaboration), Nucl. Phys. A A757, 102 (2005); K. Adcox, et al. (PHENIX Collabo- 
ration), Nucl. Phys. A A757, 184 (2005); 

[2] D. Molnar, M. Gyulassy, Nucl. Phys. A697, 495 (2002), erratum - ibid A703, 893 
(2002). 

[3] P. Huovinen, PF. Kolb, U.W. Heinz, PV. Ruuskanen, S.A. Voloshin, Phys. Lett B503, 
58 (2001); P. Huovinen in Quark-Gluon Plasma 3 eds. R.C. Hwa and X.N. Wang, World 
Scientific, Singapore (2004). 

[4] D.T. Son, A.O. Starinets, Ann.Rev.Nucl.Part.Sci. 57, 95-118 (2007). 

[5] P. Kovtun and L. G. Yaffe, Phys. Rev. D 68, 025007 (2003). 

[6] J. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998). 

[7] E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998). 

[8] S.S. Gubser, I.R. Klebanov, A.M. Polyakov, Phys. Lett. B 428, 105 (1998) 

[9] G. Policastro, D. T. Son and A. O. Starinets, Phys. Rev. Lett. 87, 081601 (2001). 

[10] D. T. Son and A. O. Starinets, JHEP 09, 042 (2002). 

[11] G. Policastro, D. T. Son and A. O. Starinets, JHEP 09, 043 (2002). 

[12] G. Policastro, D. T. Son and A. O. Starinets, JHEP 12, 054 (2002). 

[13] C. P. Herzog, JHEP 0212, 26 (2002). 

[14] C. P. Herzog, Phys. Rev. D 68, 024013 (2003). 

[15] P.K. Kovtun and A. O. Starinets, Phys. Rev. D 72, 086009 (2005). 

[16] J. Mas and J. Tarrio, JHEP 0705, 036 (2007). 

[17] P. Kovtun, D. T. Son and A. O. Starinets, JHEP 10, 064 (2003). 

[18] O. Saremi, |arXiv:hep-th/0703170 

[19] J.I. Kapusta, T. Springer, Phys. Rev. D 78, 066017 (2008) 

[20] M. Natsuume and T. Okamura, Phys. Rev. D 77, 066014 (2008). 

[21] N. Iqbal, H. Liu. larXiv:0809.3808l [hep-th] 



12 



[22 
[23 
[24 
[25 
[26 
[27; 

[28 
[29 
[30 
[31 



A. O. Starinets, larXiv:0806.3"797l [hep-th] 

D. Kharzeev, K. Tuchin, JHEP 0809, 093, (2008). 

F. Karsch, D. Kharzeev, K. Tuchin, Phys. Lett. B 663, 217 ,(2008) 

S. S. Gubser, A. Nellore, S. S. Pufu, F.D. Rocha, larXiv: 0804 .19501 [hep-th] 

S. S. Gubser, S. S. Pufu, F. D. Rocha, JHEP 0808, 085 (2008). 

S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General 
Theory of Relativity, Wiley & Sons, New York, 1972. 

A. Buchel and J. T. Liu, Phys. Rev. Lett. 93, 090602 (2004). 

A. Buchel, Phys. Lett. B 663, 286 (2008). 

B. Batell, T. Gherghetta, Phys.Rev.D 78, 026002 (2008). 

W. de Paula, T. Frederico, H. Forkel, M. Beyer, larXiv:0806.3 830 [hep-ph] 



A Background Ricci Tensor 

An explicit computation of the Ricci tensor for our background gives 

F (r) = g M R O0 = -^=d r U/=gg rr V L \g 00 ]), (67) 
F x (r) = g u R vl = —t=d r (V=gg rr V L lg xx ]), (68) 



F r (r) = g rr R rr = -L<9 r (g m g rr V L [g m ] 2 ) + /r^ (g* x g rr T> L [g xx } 2 ) ■ (69) 
4 #oo v ' W xx v 1 

Here i denotes the spatial coordinates X\...x v , and we are using the notation V L to denote 
the logarithmic derivative as defined in the text (jl]). 

Combinations which appear frequently in the text are 

F - f x = -Ld r (V=gg rr v L [f}) , (70) 



and 

fx/gTr 



F -F r = ^g rr V L [g xx ]V L 



v L [g, 



(71) 



where / is defined as in the text ([6]) 
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B Linearized Perturbation Equations 



Hydrodynamic behavior of the dual gauge theory is accessed by introducing perturbations 
on top of the black brane background. Turning on the sound mode perturbations fl24H28p . 
and expanding the background equations of motion to first order in the perturbation results 
in the following set of equations. We use the superscript "1" to denote quantities that are 
O(V)orCW). 



G 



(i) 



00 



-Jdr 



-99 



f 



-((p-l)B' + C) 



-87rGp +2 T o ) 
dU 



+ 9rr-7-m + 0(6<f>)' - (p - l)q 2 g rr g xx B = (72) 



p—i p— l 



fV=g 



zd r 



^gg rr fA> 



i=l 



-9 



dU 



-gg rr ((p - 2)B' + C) + g„ — (50) + 0'(50)' (73) 



-g rr (q 2 g xx (A + (p-2)B)+w 2 g 00 ((p-2)B + C)+2qwg xx D) = 



G« = -8ttG p+2 T« 



fy/=9 



Or 



9rr 



dZ7 



-^"■(p-l)S'J (<ty) + 0W)' 

-(p-l^W*^ = 0(74) 



^(i) 

^Oz 



-#Gr) p+1 



P 



-d r 



(g 



-g v 



-87rG p+ 2ToP 



2^ r (F - + gw(p - = (75) 



GW = -87rG p+2 T« 



dU 



V L [(g xx y}A' + V L gooig^y- 1 ((p - 1)B' + C) - 2(f>'(8<f>y + 2g rr ^-(8<f>) 
1 J dcp 

-2g„ (w 2 g 00 ((p - l)B + C) + q 2 g xx (A + (p - 1)5) + 2^^D 



(76) 



r (i) 



r 0r 



U>/<9 r 



/ 



((p-l)fl + <7) 



irrGp^T^J 



qf 2 D' + w<f/(6<f>) = 



(77) 
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w 



f 



P 



?-d r if A] + —d r f 2 D + q(p - l)B' + = 



(7f 



□ 



-5 



d 2 U" 



where we have defined 



H(r) = A{r) + (p - l)B(r) + C(r). 



(79) 
(80) 



C Derivation of Gauge Invariant Equations 



Here we present the explicit combination of the Einstein equations that lead to the gauge 
invariant equations. Consider the following combination of the Einstein equations (1721 - 179]) : 



P 



a 



P 



(JT2D - dTj) + (J7JD + iflTSDl + [dZ2D + IflTJ) 



-Iqw^h 



Here a is defined as in the text (I55j) . After a long calculation, one can show that this reduces 

to 

Z x - 2g rr {F - F x ) (w 2 {B - C) + Z ) - Bq 2 f 2 A 1 = 0, (82) 
where Z x is the left side of ( 1511) . and 



Ai = 2X? L [ fc ]a 2 <9, 



-(F -F x ) 



4#r 



;(F -F x )((f/) 



t\2 



(83) 



Clearly, Ai vanishes by the background equations of motion ( 1211) and ( fl9l) . leaving a differ- 
ential equation involving only the gauge invariant variables, namely, AMI) . 
To derive the other gauge invariant equation, the relevant combination is: 



1 + 



(Ei - (J7JD - dTJ) - i(J7j) + -<9 r 
2 J a 



w 
P 



dZZD - gdZHD . (84) 



Another lengthy calculation reduces this combination to 
2 



Zo 



V) 2 -2g rr (F -F r ))V L 
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- 5A 2 = 0, 



15) 



where Z 2 is the left side of ( 135|) . and 



A 5 



+ 



+ 




16) 



Again, all terms except for Z 2 vanish by the background equations of motion, leaving us 
with (1551). 

Generalization of these equations to multiple scalar fields is quite straightforward. In 
this case, the Einstein equations (J72] - [79]) are modified, but the combinations that lead to 
the gauge invariant equations remain unchanged. As mentioned in the text, there will be an 
additional gauge invariant equation for each additional scalar field. In this case, there will 
be n combinations like (1841) . Each will have (1T9"j) replaced by the analogous equation for each 
particular scalar field. 



D Special Case Solution 



In this section, we present the remainder of the calculation which leads to equations ff49|) . 
One must go back to (134"1) . insert = and the incoming wave ansatz (jH]), and expand the 
resulting equation in powers of q. In doing so, one should take into account the constraints 
fT2Tj) and (jUJ). The first of these constraints allows one to eliminate yf^gg rr in favor of 
T^l[P}- Once this is done, the lowest order equation for Y can be written 



(go - P^y 

a 2 f 2 



Or 



1 



[V L [P] (C -/ a a a ) 
Co 



-Y' 

2 2 



v L [P 



Co - Pa 2 



Y n = 0, 



17) 



where 

pw\-a Q . 

The general solution to that equation contains two arbitrary constants k\ and k 2 , and can 
be written 



Yo(r) = (pwl ~a + f{r) 2 a 2 ) 



k + k\ 



pw\- ap- a 2 f(r') 2 \ f'(r') 
pwl - a + a 2 f{r') 2 ) f(r') 



dr' 



19) 



as can be found by first guessing a solution Y (r) = b + bif(r) 2 , and then using the technique 
of reduction of order once this solution is found. 

The integral in (|89p is logarithmically divergent near the horizon, and thus the assumption 
of regularity leads to k\ = 0. Finally, applying the Dirichlet boundary condition at r —>■ oo, 
leads to 



pw 1 - a + a 2 = 0, 



(90) 



16 



where we have assumed that f(r — > oo) = 1. Thus, we find the lowest order term in (jUJ). 

Proceeding now to the next order in q, and substituting the solutions for w\ and Y , one 
finds the following differential equation for Y\. 



Or 



2\2 



+ 



2/f 



1+/ 



2^3 



Yi - wik 



%a 2 



+ 2pw 2 



0. 



(91) 



The solution to the homogeneous part can be found using the same techniques as for the 
solution for Y listed above. A particular solution to the inhomogeneous equation is obviously 
Yi(r) = Constant. This leads to the general solution, 



Yi(r) = (f(rf 



(/(r ,)2 _ 1)2 /(r /) 



( ia 2 

wr 



The integral above, is again logarithmically divergent, leading to the requirement that k% 
Applying now the Dirichlet boundary condition at r — > oo leads to 



(ia 2 

wr 



+ 2pw 2 



0. 



(92) 
= 0. 

(93) 



Equation (|48|) immediately follows from (|90|) and (1931) . 
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